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ABSTRACT 

Hubble Space Telescope (HST) Wide Field Planetary Camera 2 ( WFPC2) images of the Virgo Cluster dwarf elliptical 
galaxy VCC 128 reveal an apparently double nucleus. The two components, which are separated by ~ 32 pc in 
projection, have the same magnitude and color. We present a spectrum of this double nucleus and show that it is 
inconsistent with one or both compon ents being emission -line background objects or foreground stars. The most likely 
interpretation is that, as suggested bv lLauer et al.l ( 1996) for the double nucleus in NGC 4486B, we are seeing a nuclear 
disk surrounding a supermassive black hole. This is only the second time an early- type dwarf (dE/dSph) galaxy has 
been suggested to host a SMBH. 

Subject headings: galaxies: dwarf — galaxies: individual (VCC 128) — galaxies: nuclei — galaxies: 
photometry — galaxies: stellar content 



1. INTRODUCTION 

The centers of galaxies often host supermassive 
black holes (SMBHs) with masses M» ranging from 
~ 10 6 to - 10 9 M© <|Kormendv fc Richstond Tl995j) . 
SMBH masses exhibit a variet y of scaling relations 
including the M, — a re lation ijGebhardt et al.l l200Gt 
IFerrarese fc Merrittl 120001) and an M. - M n M relation 
where Mnu is th e mass of the host dark halo l)Ferraresd 
12001 IB aes etai]l200^ . These correlations are evidence 
for a coupling between the formation of SMBHs and 
galaxies. One key to unraveling this coupling is to iden- 
tify the seeds of SMBHs in dwarf galaxies. Ground- 
based stellar kinematic d ata of dwarf elliptical galaxies 
(dEs) llCeha et a,lJl2fiol rule out M. > 10 7 M . In 
the M31 globular cluster Gl, thought to be t h e strip ped 
nucleus of an accreted dE, iGebhardt et al.l l|2002|) re- 
port finding an intermediate mass black hole (IMBH) 
with M. ~ 10 4 M Q , consistent with the extrapola- 
tion of the M. — a relation. Nuclear activity in low- 
luminosity dwarf galaxies has turned up a number of 
IMBHs with estimated masses from 10 4 to 10 6 M^ 
llFilippenko fc Holl200& iBarth et aHl2004i iGreene fc Hoi 
|2004|) . The re is instead no ev i dence of a SMBH hi 
either M33 iMerritt et all 120011: IGebhardt et all 12001(1 
or NGC 205 llValliin et a.U 120051). The pa ucity of 
SMBHs in dwarf galaxies led IFerrarese! l|2002|) to sug- 
gest that low mass galaxi es are inefficient at forming 
SMBHs. iCote et al.l (12006^ . IFerrarese et"aT1 (f2006h and 
iWehner fc Harris! l)2006l) find instead a continuity in scal- 
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ing relations, with M. replaced by McmOj the mass of 
a central massive object, which is either a SMBH or a 
compact nucleus. 

HST does not permit many dynamical M. measure- 
ments below ~ 10 6 Mq. An alternative approach to 
finding SMBHs in low surface brightness dwarf galaxies 
is a morphological one. In a few galaxies, HST has found 
double nuclei; the clearest examples are in M31 and in 
NGC 4486B l|Lauer et al.l 119931 Il9961 I2005D . The two 
nuclei in M31, dubbed PI and P2, are separated by 0'.'49 
or 3.6 pc and have unequal surface brightness, with P2 
having the lower surface brightness but sitting very close 
to the global photocenter. In contrast, in NGC 4486B 
the two nuclei are equal in magnitude, color and dis- 
placement from the photocenter, with a total separation 
of 12 pc. The generally accepted model of the nucleus in 
M31 identifies P2 as the center of the galaxy surround- 
ing its SMBH and PI as the br ight off-center a poapsis 
of an eccentric Keplerian disk ( Trcmaind 119951) . Self- 
consistent mo dels and simulations of such nuclei can be 
constructed llSalow fc Staffer! 20011 Dacobs fc Sellwoodl 



l2T)0^ISamb"hus fc Sridharl2002t) . ILauer et all 1119961) oro- 
posed a similar model for the double nucleus in NGC 
4486B with the difference that the SMBH is located on 
the global photocenter between the two peaks. They sug- 
gested that binary nuclei are morphological indicators of 
SMBHs. 

Here we present another example of a galaxy with a 
double nucleus, the Virgo Cluster dE VCC 128. Its nu- 
cleus is very similar to that of NGC 4486B. Section [21 
presents the archival HST data, Section [3] presents new 
spectra and their modeling and in Section f|] we propose 
that the double nucleus is best explained by a disk sur- 
rounding a SMBH at the center of VCC 128 and explore 
its likely position in the M, — a plane. 

2. ANALYSIS OF HST ARCHIVAL IMAGES 

Using the HST archive, we explored the nuclear mor- 
phologies of the sample of dEs imaged in the WFPC2 
snapshot surveys Prop-ID GO-8600 and GO-6352 (PI 
Ferguson). These surveys imag ed 50 galaxies, which were 
analysed by lLotz et aU (|2004j) to show that typical dE 
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stellar envelopes are 0.1-0.2 mag. redder in V — I than 
their nuclei. In most of the galaxies we either found no 
nucleus or a single nucleus. VCC 1107, FCC 208 and 
VCC 128 contained what appeared to be double nuclei; 
of these the two nuclei in VCC 128 are the closest, and are 
most similar in magnitudes making this galaxy a prime 
candidate for spectroscopic follow-up. 

VCC 128 (UGCA 275) is a dE galaxy of m B = 15.6 
(Mb = —15.5) in the outskirts of the Virgo Cluster. 
VCC 128 was imaged for 460 seconds in the F555W (V- 
band) filter and 300 seconds in the F814W (J-band) fil- 
ter. Its two nuclei are separated by 4 pixels of WF3 (see 
Figure P), i.e. ~ 0'.'4 or 32 pc (assuming a Virgo Clus- 
ter d istance of 16.5 Mpc; iTonrv et alll2001t Derien et al.l 
2004). We label the component to the south-west as 
PI and the one to the north-east as P2. These are 
equally bright, with apparent magnitudes ~ 22.56 ±0.05 
(M v = -8.5 in PI) and ~ 22.72±0.06 (P2) in the F555W 
filter and V - I = 0.9 ± 0.07 (PI) and 1.15 ± 0.07 (P2). 
The S/N = 4.3 - 5.0 in V and 6.2 - 6.6 in / for the two 
nuclei measured with an aperture of radius 2 pixels. The 
V — I color map (Figure shows that the double nu- 
cleus is not caused by patchy obscuration. Each nucleus 
is resolved, with FWHM (from a Gaussian fit) of 3.1 and 
3.9 pixels in the /-band, whereas the filter point spread 
function FWHM is 1.3 pixels. 

The surface brightness profile of VCC 128 is quite flat 
and e llipse fitting in I RAF is unstable. Using GAL- 
FIT l|Peng et alJl2002f) we fitted the surface brightness 
by a Sersic model, after masking out stars and the nu- 
clei. We obtained a Sersic index n — 0.35 and a half 
light radius R eff = 9" for F814W and n = 0.55 and 
R e ff = 14'/5 for F555W. The circum-nuclear region 
(measured within 1" of the nuclei) has a surface bright- 
ness \xy — 21.37 ± 0.07 and color fiy — Hi — 1-0 ± 0.1. 
The average surface brightness within R e j / (averaged be- 
tween V and /-bands) is [iy — 23.95 ± 0.02 mag/arcsec 2 
and hi = 22.94 ± 0.02 mag/arcsec 2 . The midpoint be- 
tween the centroids of the two nuclei (measured within 
an aperture of 2 pixels) is offset by ~ 0'.'8 (~ 60 pc) from 
the photocenter of the galaxy as fitted by GALFIT but 
with these shallow images we cannot determine whether 
the photocenter varies with radius as in other offset nu- 
clei (iDe Rijck e fc De battistall2004j) . Possibly because of 
this offset. lLotz et alJ i|2004T) classified VCC 128 as non- 
nucleated. 

3. SPECTRAL ENERGY DISTRIBUTION 

3.1. Observations and Reductions 

We obtained long-slit spectra through the nucleus of 
VCC 128 at the Apache Point Observatory (APO) 3.5- 
m telescope using the Double Imaging Spectrograph with 
0'.'4 pixels. We used the low resolution blue grating with a 
resolution 2.4 A/pixel over the wavelength range 3660 < 
A < 5600 A. The data were obtained during two runs. In 
the first run on March 8, 2005 we obtained a total of 5.3 
hours of data using a l'/5 slit. Seeing conditions varied 
significantly over the night, ranging from 0'.'8 to > 2" and 
we had difficulty pointing at the nucleus during some of 
the exposures. Unsurprisingly, the double nucleus can 
be detected only in 2 one-hour exposures from this run. 
Therefore we also obtained data during a second run on 
April 21, 2006 under sub-arcsecond seeing. During this 



run we used the 0'.'9 slit to obtain 2 one-hour spectra. 

The spectra were reduced with IRAF using stan- 
dard tasks. Flux calibration, using the standard stars 
BD182647 for the first run, and BD262606, Feige 98 and 
HR 4554 for the second run, achieved a 7% uncertainty at 
wavelengths bluer than 4000 A, increasing to 10% at 5500 
A. Data from the two runs were combined by smoothing 
the 2006 data to match the FWHM of the 2005 run (6A, 
or 454 kms^ 1 at the Call lines). The spectrum of the nu- 
cleus (spatially unresolved from the ground) is typically 
6 pixels wide along the cross-dispersed axis, and that of 
the galaxy 51 pixels. We extracted the spectrum of the 
nucleus from the 4 exposures, integrated over the 6 cen- 
tral pixels. The galaxy's spectrum was extracted below 
and above the nucleus for a total extension of 45 pixels 
along the spatial direction. To remove the galaxy con- 
tamination in the spectrum of the nucleus, we subtracted 
the scaled galaxy spectrum from it. The individual spec- 
tra of the nucleus and the galaxy were then averaged 
to increase their S/N ratio, which turned out to be ~ 2 
and ~ 8 per pixel respectively. The spectra of the galaxy 
and of the nucleus are shown in FigureHJ At the nucleus, 
~ 12% of the light comes from the nuclei, the rest being 
due to the galaxy. 

3.2. Spectral analysis 

The absence of emission lines in the nuclear spectrum 
rules out, at high confidence, that either PI or P2 is a 
background emission line object. A foreground star can 
be ruled out on several counts. The nuclei are both re- 
solved by HST and stan dard models of the Milky Way 
l)Bahcall fc Soneiralll980[) predict only 0.22 stars of this 
magnitude to fall within 2R e ffy = 29" in this part of 
the sky. We also used our spectra to measure the veloc- 
ities of the galaxy and of the nuclei from the centroid of 
a Gaussian fit to the Call H and K lines. For the galaxy 
we obtained V — 1331 ± 184 km s -1 and for the nuclei 
1179±185kms~ 1 (heliocentric). As these lines have very 
low S/N, we also measured velocities by modeling the 
entire spectral energy distrib ution (SEP) as suming vari- 
ous stellar templates from the lPicklea l)1998l) library. We 
considered templates of solar metallicity, although fits 
to non-solar metallicities will not change the outcome of 
the analysis. Figure shows the reduced x 2 f° r 704 de- 
grees of freedom of a comparison in the spectral range 
3880 < A < 5000A, which includes the 4000A break, 
the G-band step, Balmer features and metallic absorp- 
tion lines. Wavelengths outside this range have very low 
S/N. Galactic reddening was appl ied using the presc rip- 
tion (R = A V /E(B - V) = 3.1) of lFitzpatrickl (tl99 9) for 
a color excess of E(B-V) = 0.015 ijSchlegel et alJ l998). 
The dots in Figure span the temperature range from 
05 (right) to M6 for dwarfs (black lines) and from 08 to 
M9 for giants (gray line) . The x 2 is shown as a function 
of distance, which is computed by comparing the appar- 
ent and absolute magnitudes corresponding to each star. 
Overall, the best fitting foreground star would be a G8V 
dwarf at 25 kpc (x 2 = 1.23) with V = 1190±f|^ kms" 1 . 
All fits with x 2 < 2 require v > 500kms _1 , which is im- 
probable for stars in the Milky Way halo. Hence, neither 
source at the center of VCC 128 can be a foreground star 
but in good agreement with Virgo Cluster membership, 
leading us to conclude that this is the nucleus of VCC 
128. 
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Fig. 1.— The WFPC2 F814W image of VCC 128 (left panel, 32" X 37"). In the center top panel we zoom into the central 2'.'5 X 1'.'5, 
showing that the nucleus is resolved into two components. PI is the nucleus to the south-west (left) and P2 the one to the north-east (right). 
The nuclei are located at RA 12 h 14 m 59 s .71 Dec 09°33'55('57 (PI) and RA 12' i 14 m 59 s .73 Dec 09°33'55'.'82 (P2) in J2000.0 coordinates. 
The center bottom panel shows the same zoom into the V — I color map, where (V — I) ranges between 1 and 2 mag, reaching ~ 1.0 in the 
double nucleus. The top right panel shows the light profile of a slit 5 pixels across along the major-axis of the nucleus. The shaded region 
indicates ±l<x error. The bottom right panel shows the color profile for the same slit. 

To estimate the mass of the nucleus, we compared the 
SEDs of both the galaxy an d the nucleus with the pop - 
ulation synthesis models of iBruzual &: CharlotJ {2003). 
We generated a 64 x 64 grid of r- models, i.e. compos- 
ite stellar populations according to an exponentially de- 
caying star formation rate that are then modeled at the 
resolution and pixel sampling of the observations, and 
compared via a maximum likelihood analysis. The for- 
mation epoch was fixed at zp — 3. The two parameters 
explored in the grid are the exponential star formation 
timescale (—1 < log(r/Gyr) < 1) and the metallicity 
(-2 < log(Z/Z©) < +0.3). Figure □ shows the 1, 2 and 
3cr (thick lines) confidence levels for the analysis of the 
galaxy (left) and of the nucleus (right) . The plots overlay 
contours of stellar mass corresp onding to the ob served 
apparent magnitude assuming a Chabrieil pOOl initial 
mass function (IMF) , and using th e total apparent mag- 
nitude of the galaxy (ros = 15.6 Bing geli et al.lll985|) 
and for the double nucleus combined. Our results for the 
ages and metallicities are independent of the assumed 
IMF but the masses are more sensitive to this change. 
The nucleus is consistent with having a relatively metal 
rich stellar population older than 8 Gyr. With the as- 
sumed IMF, the best-fit combined nuclear stellar mass 
is ~ 10 6 Mq. The galaxy spectrum gives - within error 
bars - a similar value for the age and metallicity. 

In order to assess the effect of dust on the SED analysis, 
we ran a few sets of grids for various values of E(B — 
V). The best fits obtained for these realizations gave 
consistently higher values of x 2 , so that a color excess 
above E(B — V) > 0.1 mag. is ruled out at the 3a level. 

4. DISCUSSION AND CONCLUSIONS 

Since the nucleus is old, the present configuration could 
not have resulted from recent gas infall. It is unlikely, 
though not impossible, that the two nuclei are both glob- 
ular clusters 32 pc apart in the last stages of merging 
because of the short lifet ime of such a configura tion (~ 5 
Myr using the standard iChandrasekharl l) 19431) formula 
with In A = 4 and v c = 20kms~ 1 ). Two well-separated 
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Fig. 2. — Top: the observed (smoothed) SEDs. The two panels 
compare the SED of the galaxy (left) and the nucleus (right) (gray 
line), with the best fitting model SED (black lines). The feature at 
~ 4800 A in the observed spectra is due to an artifact of the flux 
calibration due to a 'bump' in the grating response, while that at 
~ 4600 A in the nuclear spectrum is a remnant of the background 
subtraction. Bottom: the corresponding mean properties of the 
best-fitting stellar populations. Confidence levels (black lines) are 
given at the 1, 2 and 3c (thick lines) levels. Contours of stellar 
mass are shown in both cases (gray lines). For the nucleus, the 
mass shown is for the two components combined. 

globular clusters projected so close to each other, so close 
to the center of the galaxy, and so similar in their prop- 
erties is even more unlikely. The double nucleus in VCC 
128 is similar in many ways to that in NGC 4486B. As 
in NGC 4486B, a nuclear disk surrounding a SMBH pro- 
vides the best explanation: this would be a stable, long- 
lived configuration, and would account for the similar 
colors and magnitudes of the two components and for 
their location near the center of the galaxy. 

If the double nucleus is a disk orbiting a SMBH, we 
can estimate a lower limit to M, by assuming that it is 
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Fig. 3. — Top: Analysis of the nuclear spectrum assuming it is 
a star. The black lines compare the observed SED with a range of 
dwarf stars (O-M). The gray line corresponds to a range of giant 
stars. We use V = 22.6 to derive distances. The dashed and dot- 
ted lines correspond to an assumed velocity = and 500 kms , 
respectively, whereas the solid lines are computed for the best es- 
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1190 j^gQ kms 1 at 95% confidence level - computed 
by cross-correlating the observed SED with the stellar template 
that gives the best fit (G8V, M v = +5.6). The best fit for a giant 
is a G5III (M v = +0.4) with \ 2 = 1.52. Bottom: The estimated 
location of VCC 128 (star) on the M. — a plane . The solid cir- 
cles are the compilation of Trcmainc ct al. (2002), with the solid 
line showing their best fit relatio n. The dashed and dotted lines 
show the rela tion as derived by Merritt & Ferrarese] 1200111 and 
IWvithd I2006B. respectively The open circle shows NGC 4486B 
from Kormcndv & Gcbhardt (2001) and Kormendy (priv. comm.). 

larger than the nuclear disk mass, = 1O 6 M . Al- 
ternatively we ca n use the ratio of M^, to M. in M31 
(Md/M. = 0.16 iTremaind Il995l) and in NGC 4486B 
(IvW M. = 0.019 iKormendv et all Il997t lLauer et all 
1996) to estimate M.-6xl0 6 - 5 x 10 7 M Q in VCC 
128. We estimate a velocity dispersion, a, for VCC 128 
from the Faber- Jackson relation of dE's (|de Riicke et al.l 



|2005|) . We obtain a ~ 35 - 65kms _1 . Thus within large 
uncertainties the postulated SMBH in VCC 128 could 
s atisfy the M, — a relation, as we show in Figure [31 

Fer raresd l)2002fl proposed that galaxies with a dark 
matter virial velocity v V i r < 200 km s -1 would not be 
able to form SMBH's. Using her scaling relations be- 
tween a and v V i r , we find that zw < 65 kms -1 assum- 
ing the a estimated above. It is unlikely that there is 
enough scatter in the scaling relations as to accommodate 
v vir ~ 200 kms" 1 in VCC 128. Thus its SMBH proba- 
bly violates the proposed limit unless VCC 128 has been 
heavily stripped in the cluster environment. Curiously 
however, VCC 128 would still satisfy th e relation between 
Mcmci and -B— band magn itude Mb <|Wehner fc Harris! 
2006: Jtorese et al.ll2006|) . 

We have shown that the nucleus of VCC 128 is double, 
with the two components of equal magnitude and color. 
We proposed that this can be explained by the presence 
of a disk surrounding a SMBH. Simple estimates of its 
mass all put it in the regime M, ~ 10 6 M Q or larger. 
If VCC 128 is on the fundamental plane, then these es- 
timates put the SMBH at, or above, the M. — a rela- 
tion. However the halo virial velocity would be smaller 
than the proposed limit for galaxies to be able to form 
SMBHs. Amongst the dwarf elliptical/spheroidal pop- 
ulation of galaxies, this is only the s econd example of 
a sys tem with evidence for a SMBH ijMaccarone et all 
2005|. Given the very interesting nature of this object, 
we suggest that further high-resolution imaging and spec- 
troscopy of VCC128 would be well worthwhile. 
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